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SUMMARY 



The development of the laminar TDOundary layer in a 
comprasslTjle fluid is considered,- Formulas are given 
for determining the •'boundary-layer thickness and the ratio 
of the houndary-layer Reynolds numher to the hody Reynolds 
number for airfoils and "bodies of revolution. 

It is shown that the effect of compressibility will 
profoundly alter the Reynolds number corresponding to 
the upper limit of the range of the low^-drag coefficients. 
The available data indicate that for low-drag and high 
critical compressibility speed airfoils and bodies of rev- 
olution, this effect is favorable. 



INTRODUCTION 



Experiments with a large number of low-drag airfoils 
have shown that as long as- the transition from laminar to 
turbulent flow fit the surface occurs beti^een the minimum 
pressure position and the trailing edge of the airfoil, 
the -low-drag Characteristics of these airfoils aro main- 
tained, but that as the transition point moves forward of 
tho minimum pressure position the drag coefficient in- 
croa.sos more or loss markedly depending on tho airfoil 
pressure distribution. 

It hp^e been found thnt the boundary-layer Reynolds 
number Rg , based on the boundary-layer thickness and 
the local velocity outside this layer, gives a f^ir meas- 
ure of the stability of the boundary layer and, in conse- 
quence, may be used n s a criterion for determining the 
point at which transition to turbulent flow takes place. 



As noted in reference 1, the best estimates of the 



crit icnl value of 



availp/ble present were olstained 



from flight tests of an lAGA 3d-215 nirfoil section which 
was tested as a glove on the l^-iS airplane. Designating 



V 



where V is the velocity outside the Isoundary layer; 6, 
the distance from the surface of the airfoil to a point 
in the "boundary layer v/here the velocity has reached 
0.707 T ; and U, the kinematic, viscosity of the fluid, 
critical values of Eg "between 8000 and 9500 were oh- 
served in these flight tests. 

Irom von Zarmpn's momentum relation (reference 2, p. 
107) it is evident that, in an incompressible ' fluid., if 

1. the "boundary layer on a given hody is laminar 
from the stagnation point to any given point on the "body 

2, the houndary-layer velocit;- profile is at- all 
points along the • surface of the same form when consider-Sd 
nondimen si onally in terms of 6 and - V; then at -the 
given point on the surfac-e , the 'boundary-lajrer Eeynolds 
numher Eg is related, to the "body "Eeynolds number ■ E 

("based on body dimensions and the stream speed) in the form 



E, 



E 



constant 



If the ,c.ons.t.ant is .known for any body at the minimum 

pressure point i.t i s .posiS'i.ble to determine t-he body. 

Eeynolds number, vrhich is the upper limit of the range of 

the low-,drag coefficients, -for -a given value of Ee- . . 

°cr 1 1 

For nearly inc'ompr e s sible flo\«r this constant may be. eval- 
uated' by the metliod >f reference 1. In those applica- 
tions "where the HaCh number is not negligi.bly .small",; i.t .is 
necessary to extend this method to take account of the 
compressibility effects.. Such an extension of this meth- 
od is the subject of this paper. 



TH30EY 



The growth of the laminar houndary layer in a com- 
pressible fluid may be conveniently studied by von 
Karman's momentum method. To this end, consider first 
the steady-state flovi across the faces of an elemental 
parallelepiped at the surface of a two-dimensional body 
shown in figure 1. Let h, v;hich is chosen so as to be 
independent of s, be the distance from the surface of 
the body to a poiat in the boundary layer, where the fluid 
shear has become negligibly small.. 

The several contributions to the s component of the 
change in momentum across the parallelepiped *will now be 
considered in turn. 

The fluid entering the face normal to s per unit 
width introduces the momentum 




p 11 dy 



while that removed at the opposite face is 



Jph ph 
p u^ dy+ d6 / p u" 
n '-J n 



dy 



and hence the change in this contributipn to the momentum 
.is ■ 



ds 



Lis J p " 



dy 



(1) 



Ho contribution occurs at the surface of the airfoil but 
at the parallel face an amount py 7^ ds is removed. 
Continuity requires thnt 



d(pv) 
9 y 



3(pu) 
B s 



hence 



4 



^ = - J, p ^ 

and so this cont JSi"but ion becomes 

h 



- V ds 



/ p ^ 



(2) 



Since the flow is considered to be constant with time, the 
total change in the -s component of momentum across the 
parallelepiped is given by the sum of equations (l) and 
(2). ■ 

The forces acting on the parallel opipcd in the direc- 
tion of s are the stxrface shear and the pressure dif- 
ference between the surfaces normal to s. The shear 
force, using the established sign convention, is 

• - T ds (3) 

and, if the boundary layer is thin, it has been shown 
(reference 2, p. 83) that the pressure variation with y 
is negligible, so that the pressure force is 



ds 



- h — ds 



Now V is small compared to V, so that Bernoulli's 
equation for a compressible flow which is constant with 
respect to time may be written for the flow region, out- 
side- the frictiohal influence of the boundary layer 

dp dV ■ 1 d (p^Y^ ) 1 dp^ 

ds ds 2 • ds 2 . ds 

It is convenient ta rewrite this as 

ds . ds ds ds 

d (pttY^) d (p^v) 

= 1 _ + Y — 

ds ds 



Moreover, since "both.- p^' and V are independent of y, 
for reasons which vill "be evident later 



h h 

dp 1 d r , . V d r 

d^ = - h di I P-^^^ '^Ts I 

so that eq^uation (4) "becomes 



^•-'o "o ■ 



(5) 



finally, equating the change in the s component of mo- 
mentum across the parallelepiped to the s directed 
forces on the parallelepiped, the "momentum" relation for 
the two-dimensional flow of a compr e s si'ble fluid (i.e., 
with Varying density) is 

= ^ / ^P^^^ - P - ^ ^ j ^Pv^ ~ P u) dy (6) 



It has 'been ohserved in a numlier of experiments with 
conventional low-drag and high critical compressihility 
speed airfoils that the Blasius-type profile is a good 
approximation to the actual houndary-layer profile over the 
forward region of the airfoil where the pressures are fall- 
ing. An examination of the calculated houndary-laz/or pro- 
files for' a flat plate at a number of Mach numbers (refer- 
ence S) indicates that the form of the profile remains 
close to the Blasius type for subsonic flows. As a conse- 
quence it seems reasonable to assume, as is done in the 
analysis to follow, that the boundary layer over the sur- 
face of conventional low-drag and high critical speed air- 
foils will remain of the Blasius type throughout the sub- 
sonic speed range. 

It may be shoxirn that to the order of the density 

and temperature outside the boundary layer, if adiabatic 
conditions exist, are, respectively, 



.6 



and 



Y-1 



M' 



(7) 



where the subscript o denotes conditions in the free 
stream and the subscript v denotes conditions just 

outside' the houndai'y .layer at any point s along the 



airfoil, where the velocity is 7; and V = Cp/( 
ratio of specific heats. 



the 



Sincethe pressure is transmitted unchanged through 
the "boundary layer, it follows from the laxir of Boyle and 

Charles that the density at any point y within the 
"boundary layer is related to the local temperature "by 



P P. 



(3) 



purther, it is shown in reference 3 that for a flat plate 
the temperature varia.tion within the "bound-ay layer, for the 
Prandtl number equal to unity, is given by 



The Prandtl .number is denoted by 

k . 

wherein, for the fluid, 

the absolute yiscositj'- coefficient 



the specific heat at constant pressure 
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and 

k tho thermal conductivity 

por air, the Prandtl numher is less than unity (at standard 
condition Pr for air is 0»733) "but it is not expected 
that the form of the temperature variation as given in 
reference 3 will, for air, "be seriously in error. For the 
airfoil it seems correct then to assume the "temperature 
variation to he of the same form 



T = + { 1 - f 1 - r^Y"^ (9) 

Moreover, the results of tests with a circular cylinder 

(reference 4) have indicated that the surface temperature 
may he given with reasonable accuracy by 



^u=o = 



(10) 



Finally, from the relations of equations (7), (8), (9), 
and (10), to the order , it may he found that 



P = Po \ 



- 1 + CY-l) (Pr*) 



?he surface unit shear is given hy 



Experiment has shown that p ' varies as the absolute 
temperature to the 0o75 power and from equations (7) and 
(10) to the order of 



1 



(12) 



and so, to the pi-esent orior of approximation 



0.76 (Y-l) It 



Using tho density relations of equations (?) and (ll), 
the value of T given "by equation (IS), and assuming the 
Blasius variation of u/V vrith y/h in the momontuH 
cq^uation (6), it was found that, to the order of , tho 
hound ary-layer thickness 6 is given "by 



6^ = 



cv 



^1 VJ 



a , 1 7 r f 

5 , 3 -j^l + M 



-0.35 h- / (-1-) d(s/c)-0.44 M-^ / i j d(s/c) 

v;her e ■ • :. ' ■ ■ ■ ... ■ ' 

c chord of the airfoil 

the kinematic viscosity in the free stroam 

T]_ . the velocity outside the "boundary layer at the point, 

^i/c which the houndarj'' layer is "being computed 

6 tho "boundary— layer thickness, which is considered in 

this analysis to "be tho distance from the surface of 
the airfoil to a point in tho "bouuidary layer v/horc 
tho ratio of tho local velocity to the velocity out- 
side the "boundary layer is 0,707 



(14) 
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lo employ as a criterion for determining the 

stalDility of the laminar "boundary layer, accou.nt mxist te 
taken of the fact that "because of aerodynamic heating, 
the kinematic ■viscosity varies throughout the "boundary 
layer. The value of x> used in calculating the "boundary- 
layer Reynolds numljer should "be that characteristic of the 
point in the "boundary layer at which insta"bility initiates. 
Ihe theoretical analysis does not indicate the location of 
this point, "but it is clear that the characteristic vis- 
cosity v/ill lie "between that at the airfoil surface and 
that at the outside of the "boundary layer. The kinematic 
viscosity is given "by 




(15) 



Using the density and temperature relations of 
eq^uation (7), the "outside" viscosity is to the order 
of 



V o 



^ 2 



0 



1 - D.76 (Y-l) 



V. 



1 + 0.35 M 



V 




- 1 



(16) 



and using the relations of equat ions ( 11 ) and (l2) the 
"inside" viscosity is to the order of 



M 



1 +- 



3- Q-^J |l+(V-l)|(^Pr"^->0.76(^l-Pr"^j }-l+0. 76(Y-1 ) 



"4 



1 + M 



(0.65 (1-J 



0.35 



(17) 



10 



If . 



cV 



(18) 



then with the value of 6 given hy equation (l4) 



a. 17 



1/ c 



- 0.4 



(19) 



and the values of (a) and (h) are for the t\^o limiting- 
cases . 

(1) "based on the "outsicLe" viscosity a = 0. 02 ; h = — 0, 34 

(2) TDased on. the "inside'* viscosity a =0. 63-; h =— 0, 34 

• ; (20) 

Por the compr ossihle fluid "boundary layer of a hody 
of revolution, the momentum relation may ho found to ho 



T r = 
0 dsJo 



n 

■1 



s 2 
P^^. - P ^ 



• r dy 



h 



p y >- p u 

' V ^ 



r <iy (21) 



and under tho previous assumptions as to temperature and 
density variation and the shape of the houndary—layer 
velocity profile, it may ho shown that to the order of 
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d( s/L)- 0-» 44 U 



© 



cL(s/L) 



(22) 



and 



1/: 



' ' d(s/j.) 



-0.44 y ^: 



,N^S ^^ \10.17 



d( s/L) 



(23) 



v/hore for the two limit ing v is co s it ios , tlio values of (a) 
and ( Id) J arc given lay (SO) and 



^•1 
r 



the radius of the 1)0 dy at 

the radius of the hody at s where the velocity 
is ? 



Ii the length of the hody 

and the remaining synhols are as previously designated. 

To apply the eq_uations, the velocity dist r ihut i on at 
the Mach numher M. must he ascertained. ¥hen the exper- 
iniental pressure coefficient P distribution is known at 
the desired Mach number, the distribution of "V/Tq aay be 
found using Bernoulli's equation for a compressible fluid, 
For air this equation is 



^ Y 1 - [1+0.7025 M^P] 

0.2035 M3 



0.3 883 



(24) 
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Values oTstainod from this oquation are gi-Tcii in talolc I, 

In .the more u.sual case whore the pressure— coefficient 
distriljut ion is Icnown f or II = 0, -that for the desired 
i'iach numher may be calculated using Ton Karman' s equation 
(reference 5) 



/ l-H^ + 

2C1 + y 1-M*] 

Yalvxos ohtained from this equation. are given .in t-ahlo H» 

Discussion AiTD COFCLUSIOHS ■ 



An investigation of the ho^mdary— lay or thickness at a 
point 55 percent hehind the leading edge on the upper sur- 
face of an NACA 55,2—420 airfoil at several Hach numbers 
was conducted in the IS— foot wind tiinnel at the Ames Aero- 
nautical Laboratory, Using the measured pr es sur e cli str i-* 
biitions at the same Mach niimbers, the bovindary— layer 'thick- 
ness was calculated bj^ equation (14) v;hich considers effects 
of corapr essibility and aerodynaaic heating, and by the corre- 
sponding equation of reference 1 vrhich neglects these effects. 
The calc-alated variation of boundary— laj'-er thickness with 
liacii number as determined from these equations and the sev- 
eral experimentally—measured values are shown on figure 2-, 
That the theoretical variation of 5 is valid is inQ.icatod 
by the close agreement between the calculated va-luos obtained 
from equation ■■ (14) 'and. the experimental results. 

As noted previously, in order that Hg may bo used 
as a criterion for the stability of the boundary layer, it 

is essential, to determine where, vithin the boundary . lay or , 
transition to turbulent f 1 ow in it iates , In oxp or imont al 
investigations with flap pl-ates (reference 6), -it was found 
that slow fluctuations of flow occur within the boundary 
.layer though they are not apparent near the outside of the 
layer, Jones (reference 7) obtained experimental results 
substantiating these data and stiggested that the phenomena 
of transition to turbulent flow may be the direct result of 
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intermittent instability due to transient separation 
of the flow from the surface. If this is true, then 
transition must initiate near the inside of the "boundary 
layer, 

Ihg'exper imental data availalile showing the effects 
of. coapr essihility and aerodynamic heating, indicate that 
transition does arise near the surface'o These data were 
oh'tainred with four HACA 27—212 airfoil models of different 
chords hy measuring the maximum Reynolds number for which 
lovr drag v;as maintained. It vras found that the valties of 
this critical Reynolds number for the smaller chord airfoils, 
which required higher Mach nti.mbers than the larger chord 
airfoils to reach a given Reynolds number, were much higher 
than those for the larger chord airfoils. The variation 
of critical Reynolds number v/ith Hach number computed from- 
equation ( 19 ) and experimental measurements for the HACA 
27>-212 airfoils are shown on figure 3, It is soon that 
when the boundary— lay er Reynolds number is based on the 
inside viscosity and the low Hach— number experimental 
points are fitted to the theoretical cvirve, the calculated 
effect of compressibility is in agreement with experiment. 

To investigate the Hach— number effect for other types 
of airfoils, values of Rg^/R^ have been calculated for 
the foil owing! 

1, the HAGA 35—215 section, v;hich is representative 
of the more usual low— drag airfoils; 

2, the i-TACA 16—212, which is representative of the 
thin high critical compressibility speed airfoils; and 

3, the NACA 45—125 i^uiich is representative of the 
thick high critical compressibility speed airfoils, d.eveloped 
to permit power plant installations within the vring, Por the 
first two airfoil sections the theoretical velocity distr'i-*- 
bution was employed, Jor the third airfoil section,, the 
velocity distribution obtained from the experimental pressure 
distribution given in reference 8 vras u.sed. In figure 4, the 
ratio of Rg^/R^ as calculated by oqua.tion (19) of this re- 
port, using the limiting values of. a and b given by 
eouation (20), is shown as a function of the Hach number for 
each of the three airfoils. lov comparison tho ratio as 
calculated by tho method of reference 1, which disregards 
the compressibility effects, is shown. 
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Eecently. an experimental investigation was made of .' 
the effect of heating the surface of a lo\-/— drag airfoil 
"by heating elements placed v/ithin the wing. It vras found, 
when heat was so applied as to maintain the entire sxirface 
of the airfoil over which the laminar flow occurred at a 
constant temperature increment ahove the temperature of 
the amoient stream, that the boundary layer was destahilized 
so that the critical Eeynolds number was decreased, ¥ith 
the same temperature increment hut with only the surface in 
the immediate vicinity of the transition region heated, the 
dest aoilizing effect on the "boundary la,yer was even more 
marked, These resultSj in contrast to the results of the 
experiments v/ith the STAOA 27—313 airfoils previously alliided 
to, would indicate that viscosity considerations alone are 
not sufficient to explain the effect of heat on the stability 
of. tiac boundary layer unless, as in the cases of aerodynamic 
heating and no heating, the temperature gradient in the bound- 
ary layer at the surface is zero. It was considered that, 
in the investigation of the heated lovj— drag airfoil, by heat- 
ing the loading— edge section" only a temperature variation in 
the boundary layer • at transition similar to that obtained- 
in aerodynamic heating could be promoted So that a further 
study of aerodynamic heating could be made. This was at—, 
tempted but the h eat so transferred, vras . insufficient to 
mat er ially ■ i-nf luence the -temperature var iat ion .over and 
above that occurring natxirally at the ilach numbers of the 
tests so that no further conclusions could be dravrn. 

The derivation of equation (19) is based on the assump- 
tion of small values of the Ilach number, lever the 1 as s , the 
ec[Uation seems to be in fair q\iantit.at iv e ..agreement with 
experimental data even at Hach numbers of. 0.65, Hovrever, 
no account has been taken of effects on the wind— tunnel 
data resulting from the tunnel turbulence level changing 
v/ith speed, A recent investigation in the 8— foot high— speed 
t-annel at the Langley Memorial .Aeronaut ical laboratory in- 
dicated' that the tiinnel turbulence increasi.es with the Hach 
nu.mber. Since increasing the tur bul.enc.e level has a de- 
stabilizing effect on the boundary layer (fig. 3), it would 
be concluded that the variation predicted by equation ( 19 ) 
is' conservative, further investigation of the effects of 
coaprossibility on the stability of the laminar boundary 
layer is nc-cessary before an evaluation of th.e quantitative 
accuracy of equat io-n- ( 19 ) can . be. made,,- . 

Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Hoffott field, Calif, 
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lABLS I 



Values of - yi"^ 
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0 
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,0011 
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, .0013 


,0050 
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TABLE II 
Values of (P~I^_o) 



\. M 






u,3 


n k 




0 fi 


n 7 




+0.7 


0.0022 


0.0091}- 


0,0220 


0.01J28 


0.0660 


0.1032 


0.1550 
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+.6 


.0021 


.0086 


.0201 


.0392 


.0619 


.097^ 


.1^75 


.226 
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figure 1.- Botaadary layer coordinatas. 
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Fig-are 2.- The effect of conprossibility on the Id omadary- layer thickness, 5. HACA S6,2-i20 • 
airfoil: syc = 0.55. ^3 
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ffigure 3.- The effect of compressibilitjr on the critical Reynolds- 
n-um"ber ratio of the lACA 21-7X2 airfoil. 
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CalcTilated "oy equation 19, compressibility considered, 

"outside viscosity" 
Calculated lay equation 19, comipressi'bility considered. 
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Pigure 4.- The offect of conipressiMlity on /iL^ for three airfoils. 



